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Opto-electronic devices utilizing graphene have already demonstrated unique capabilities, which are much more 
difficult to realize with conventional technologies. However, the requirements in terms of material quality and 
uniformity are very demanding. A major roadblock towards high-performance devices are the nanoscale variations 
of graphene properties, which strongly impact the macroscopic device behaviour. Here, we present and apply 
opto-electronic nanoscopy to measure locally both the optical and electronic properties of graphene devices. This 
is achieved by combining scanning near-field infrared nanoscopy with electrical device read-out, allowing infrared 
photocurrent mapping at length scales of tens of nanometers. We apply this technique to study the impact of 
edges and grain boundaries on spatial carrier density profiles and local thermoelectric properties. Moreover, we 
show that the technique can also be applied to encapsulated graphene/hexagonal boron nitride (h-BN) devices, 
where we observe strong charge build-up near the edges, and also address a device solution to this problem. The 
technique enables nanoscale characterization for a broad range of common graphene devices without the need 
of special device architectures or invasive graphene treatment. 


As large scale integration and wafer scale device pro¬ 
cessing capabilities of graphene have become available, 
technological implementations of electronic and opto¬ 
electronic graphene devices are within reach.At the 
same time, to achieve high device performance, any im¬ 
perfections at the nanometer or even atomic scale need 
to be minimized or even eliminated. For example, in 
large area graphene, grown by chemical vapour deposi¬ 
tion (CVD), grain boundaries are the stitching regions 
between different mono-crystalline parts of graphene 
and act as carrier scatterers, limiting the graphene 
mobility and uniformity.These nanoscale defects 
are elusive to many standard characterization tech¬ 
niques without special treatment of the graphene. 

In addition, even perfectly monocrystalline graphene is 
still highly sensitive to its environment, and on typ¬ 
ical substrates charge-density inhomogeneities (charge 
puddles)and additional doping near con¬ 
tacts, defects and edges arise, which reduce the device 
performance as well. Therefore it is important to effi¬ 
ciently probe the nanoscale opto-electronic properties of 
graphene and to understand the microscopic physical be¬ 
haviour. 

A major challenge is that many of the avail¬ 
able techniques are invasive,rely on specifically de¬ 
signed device structures,image only very small 
areas,rely on high doping of the graphene,^ 
need unhindered electrical access of the probe to the 


graphene,or lack the desired nanometer 
resolution^^ and are difficult to implement. Therefore, 
a nanoscopic tool that probes both electrical and optical 
response of graphene devices at nanometer length scales 
is highly desired. 

Here we demonstrate fully non-invasive room- 
temperature scanning near-held photocurrent 
nanoscopy^^“^^ for the hrst time applied with in¬ 
frared frequencies and use it to study the nanoscale 
opto-electronic properties of devices that can later be 
used for real applications. This technique allows measur¬ 
ing the properties of graphene devices that affect their 
performance with high spatial resolution in atmospheric 
conditions. We apply this technique to study the micro¬ 
scopic physics of grain boundaries and charge density 
inhomogeneities. In addition, we study encapsulated 
graphene devices,where the encapsulation would 
prevent many other scanning probe techniques from 
accessing local properties of graphene. In general, this 
technique operates most effectively with mid-infrared 
light because it does not lead to photodoping^ and it is 
more stable in operation, compared to visible light. 

The measurement principle is sketched in Fig. la. The 
setup is based on a scattering-type scanning near-held 
optical microscope (s-SNOM)^’^^ augmented with elec¬ 
trical contact to the sample to measure currents in situ. 
In contrast to conventional s-SNOM we do not need to 
measure the outscattered light but rather directly mea- 
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Figure 1 | Near-field photocurrent working principle and photocurrent from grain boundaries, a, Sketch of the scattering- 
type scanning near-field optical microscope setup. A mid-infrared laser illuminates the atomic force microscope tip, which generates 
a locally concentrated optical field, which is absorbed by the graphene generating a position dependent photocurrent. The blue re¬ 
gion in the graphene lattice represents a grain boundary with a modified Seebeck coefficient. The arrows sketch the photocurrent 
flow pattern. For each position only the magnitude and direction of the current are measured. The sketch is not to scale, b, /pc 
map at at backgate voltage Vbg = 0 V of a CVD graphene device with three contacts: top left (drain), right (source) and bottom 
left (ground). Both grain boundaries and wrinkles show characteristic photocurrent patterns. The green box indicates the measure¬ 
ment region in c,d. c. Topography of etched CVD graphene does not show grain boundary but only wrinkles and other inhomo¬ 
geneities due to the transfer process, d, /pc at Vbg = 0 V clearly shows a grain boundary and the expected sign change around 
it. The black dashed line indicates the measurement positions in Fig. 2a, d. e. Topography (brown) and /pc (red) measured at the 
brown line in c and the red line in d respectively. 


sure current induced by the near-held as explained in 
the following. A 10.6 pm mid-infrared laser illuminates 
a metallized atomic force microscope probe, tapping at 
its mechanical resonance frequency. Part of the incoming 
light, polarized parallel to the shaft of the probe, excites 
a strong electric held at the tip apex due to an antenna 
ehect.^ The spatial extent of this near-held is on the or¬ 
der of 25 nm, limited only by the tip radius and much 
smaller than the free space wavelength of the impinging 
light. ^ 

The near and far helds impinging on the device induce 
charge hows in the device (by mechanisms discussed be¬ 
low), and drive currents into an external current ampli- 
her via contacts on the device. We isolate the part of 
the current that is induced by near helds by demodulat¬ 
ing the current at the second harmonic of the tip tapping 
frequency.^ This demodulated current is denoted Ipc and 
referred to as near-held photocurrent and is obtained to¬ 
gether with near-held optical and topography informa¬ 
tion. A typical map of Ipc , obtained by scanning the tip 
over a CVD graphene device, is shown in Fig. lb. 

We can assess the spatial resolution of the photocur¬ 


rent maps by comparing a region near the edge (Fig. Id) 
with a topographic image from the same region (Fig. Ic). 
As can be seen, Ipc falls to zero for tip locations away 
from the graphene on a similar length scale as the topog¬ 
raphy, demonstrating the successful isolation of near held 
contributions. In Fig. le we quantify the resolution by 
observing the change in Ipc as the tip is moved over the 
edge of graphene. The full width at half maximum of the 
photocurrent peak at this location is ^ 100 nm, matching 
the rise distance in the topographic signal. This resolu¬ 
tion is far below any limits relating to the 10.6 pm free 
space light wavelength. 

As to the physical mechanism of the photocur¬ 
rent, we consider the photo-thermoelectric effect that 
has been shown to dominate the photoresponse of 
graphene:^’^^^^^’^^’^^ the light (in this case, the tip- 
enhanced near-held) locally heats the graphene, and this 
heat acts via non-uniformities in Seebeck coefficient S 
to drive charge currents within the device and into the 
contacts (see Methods). Therefore, we interpret the 
variations of Ipc in terms of microscopic variations in 
S. The Seebeck coefficient, which depends on mate- 
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Figure 2 | Photocurrent profile at a grain boundary and 
its gate voltage dependence, a, Photocurrent profiie, mea¬ 
sured at the biack dashed iine in Fig. Id, perpendicuiar to the 
grain boundary at Vbg = 0 V shows good agreement with the 
photo-thermoeiectric modei with /cool = 140 nm. The modeiied 
spatiai Seebeck profiie (with FWHM 20 nm) is shown in biack. 
b, Two-probe device resistance as a function of Vbg- c, Sim- 
uiated Seebeck coefficient Sg for pristine graphene (soiid iine) 
and Sgb for poiycrystaiiine graphene with an average grain 
size of 25 nm (dashed iine) (see Suppiement, II D). d, Back- 
gate dependent photocurrent profile perpendicular to the grain 
boundary shows that the grain boundary changes its sign at 
the charge neutrality point, e, Simulated backgate dependent 
photocurrent profile based on the Seebeck profiles in c. 


rial properties such as carrier density and mobility, is 
a measure of the electromotive force driven by a tem¬ 
perature difference in a material. A complete descrip¬ 
tion of Ipc needs to take into account the carrier cool¬ 
ing length^’and overall sam ple geometry.The carrier 
cooling length /cool = \//^/^, where k the thermal con¬ 
ductivity in plane and g the thermal conductivity out 
of plane to the heat sinking substrate, describes how far 
heat propagates through the charge carriers, before dissi¬ 
pating to the environment.^^ A quantitative model of the 
thermoelectric photocurrent mechanism can be found in 
the Methods section and in the Supplement. 

We first discuss the application of this infrared near- 
held photocurrent technique to grain boundaries, which 
are responsible for some of the line-shaped features in the 
photocurrent map in Fig. lb. The region within the green 
frame is shown with higher resolution in Fig. Id, exhibit¬ 
ing a strong photocurrent signal that changes sign along 
a sharp boundary, yet the graphene is topographically 
hat in the vicinity of this boundary (Fig. Ic). We show 
now that this type of feature indicates a grain boundary. 

Figure 2a shows a line prohle of /pc across the bound¬ 
ary feature identihed in Fig. Id. This antisymmetric Ipc 
can be explained by a localized deviation in S at the 


boundary, i.e. a line defect within an otherwise uniform 
thermoelectric medium. Indeed, grain boundaries behave 
as localized lines of strongly modihed electronic proper¬ 
ties, within otherwise uniform graphene.We 
remark that the decay of the photocurrent away from 
the boundary extends over more than 100 nm, which is 
due to a larger hot carrier cooling. We hud in this case 
kooi = 140 nm. 

To gain more insight in the Seebeck coefficient at the 
grain boundary, we tune the carrier density by a global 
gate (Fig. 2d). We observe that the antisymmetric spatial 
photocurrent profile changes sign as the backgate volt¬ 
age Vbg passes the peak in resistance, i.e., the global 
charge neutrality point Vp. The Seebeck coefficient Sq 
of graphene itself changes sign at the charge neutral¬ 
ity point^’^^’^^’^^ (Fig. 2c), thus our data implies that 
the Seebeck coefficient of the grain boundary 5 'gb is al¬ 
ways smaller in magnitude than S'g, since /pg(4^bg) oc 
5'g(4^Bg) — 5'GB(kBG)- 

Using a poiycrystaiiine graphene model, we compute 
the resistance due to grain boundaries using a Kubo 
transport formalism and real space simulations.^^ 5'gb is 
the ratio of the first- and zero-order Onsager coefficients 
(see Supplement, II D). Indeed we find that S'gb is always 
smaller in magnitude and has a similar lineshape as Sq 



Figure 3 | Photocurrent from charge puddles, a, Near¬ 
field photocurrent map of an exfoliated graphene device on 
300 nm Si02 at Vbg = 20 V. The dashed lines indicate the 
position of the contacts and solid lines the graphene edges. 

The green box indicates the measurement region in b. b, De¬ 
tailed photocurrent map at the charge neutrality point of the 
device (Vbg = 7 V) reveals the charge puddles and the high 
spatial resolution of the technique, c, Autocorrelation of the 
photocurrent from charge puddles at Vd (data points) com¬ 
pared to photocurrent expected from a random charge puddle 
distribution and /cool = 200 nm (blue curve). Autocorrelation is 
taken along the source drain current path. 
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Figure 4 | Dependence of photocurrent profiles on backgate voltage reveals doping inhomogeneities. a, Backgate depen¬ 
dence of the resistance of the device measured simultaneously to the photocurrent in blue. The red curve shows the normalized 
root mean square of the photocurrent across the device. The green curve shows a single normalized photocurrent backgate trace, 
corresponding to the green dashed dotted line in c. b, Backgate dependent Seebeck coefficient of graphene, calculated from the 
gate dependent resistance in a using the Mott formula.c, Backgate dependence of the photocurrent across the device. Graphene 
is between the black dashed lines, which indicate the edges of the metal contacts. 


in the carrier density range measured (Fig. 2c). Fig. 2e 
shows a simulation of the photocurrent for the calculated 
Seebeck coefficients, which is in excellent agreement with 
the measurements. 

We next examine near-held photocurrent in a typi¬ 
cal two probe exfoliated graphene device (Fig. 3). A 
strong photocurrent is obtained with the tip near the 
metal contacts, similar to previous near- and far-held 
measurements.Additionally, an apparently ran¬ 
dom pattern of photocurrent is present throughout the 
device, as in high-resolution far-held measurements^^ but 
at a much hner scale. 

The random photocurrent pattern between the con¬ 
tacts in Fig. 3a indicates random variations in Seebeck 
coefficient over short length scales. Random variations 
of the Seebeck coefficient are indeed expected since it 
depends on carrier density,^^ which in turn has hue- 
scaled inhomogeneities (charge puddles).The 
photocurrent variations can thus be used to gain insight 
in the charge puddle distribution. A more detailed view 
of the photocurrent due to charge puddles in Fig. 3b 
shows that the length scale that can be resolved is on 
the order of hundreds of nanometers. 

Quantitatively, from the autocorrelation of the pho¬ 
tocurrent in comparison to a photo-thermoelectric model 
taking into account the size of the charge puddles in 
Fig. 3c we extract /cool ^ 200 nm. The charge puddles 
are modelled to have a size of ^ 20 nm, in accordance 
with measurements of graphene on 

By changing the gate voltage we study the carrier den¬ 
sity profile with high spatial resolution (see Fig. 4) and 
highlight the possibility of spatially resolving the charge 
neutrality point for a large device. Jpc from charge pud¬ 
dles is largest around the charge neutrality point and 
varies with position. This is consistent with the very 
high sensitivity of the Seebeck coefficient to changes in 
carrier density, near zero density (Fig. 4b). This allows 


us to map the local carrier density offset (charge inho¬ 
mogeneity) throughout the device, as indicated by the 
extremum of photocurrent in a scan of photocurrent vs. 
gate voltage (Fig. 4c). In contrast to the grain boundary 
photocurrent we do not observe the charge puddle pho¬ 
tocurrent change sign when sweeping the average carrier 
density through the charge neutrality point, as expected 
from the gate dependence the of Seebeck coefficient. 

We can thus resolve the local charge neutrality point 
at a given position of the device (green curve. Fig. 4a), 
which can be different from the global charge neutrality 
point W), the backgate voltage Vbg at which the resis¬ 
tance is maximum (blue curve. Fig. 4a). We show that 
the global charge neutrality point (blue curve. Fig. 4a) is 
determined by an average of the gate voltages at which 
the local charge neutrality points appear (red curve. 
Fig. 4a). Spatially resolved puddle photocurrent can be 
much narrower (green curve. Fig. 4a) than the average 
of all possible current paths (red curve. Fig. 4a) This in¬ 
dicates that the graphene locally has less inhomogeneity. 
Thus the technique gives insight not only in the global 
but also in the local behaviour of the device. 

Finally we apply this technique to a graphene de¬ 
vice encapsulated between two layers of h-BN, using the 
polymer-free van der Waals assembly technique^^’^^ as 
sketched in Fig. 5e. This device lies on top of an ox¬ 
idized silicon wafer, used as a backgate. The stack is 
etched into a triangle and electrically side-contacted by 
metal electrodes.Recent studies^^’^^ have shown that 
the edges affect where current flows in the device, in par¬ 
ticular near charge neutrality. In the following we study 
the build up of edge doping and provide a solution to 
this. 

While monitoring the photocurrent of such encap¬ 
sulated devices, we observe indications of strong car¬ 
rier density variations near the edges over micrometer 
scales. These variations are influenced by lighting con- 
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Figure 5 | Near-field photocurrent maps revealing edge doping in encapsulated graphene, a, Spatial photocurrent profile 
vs. backgate voltage Vbg (minus voltage of the resistance maximum Vd) near the edge of encapsulated graphene. These data are 
taken directly after annealing the device, b, The same scan on the same device after three hours in air, and c, after annealing and 
applying Vqq up to 3 V. d, The same scan after ~ 20 hours in air and after applying Vbg up to 50 V. e, Sketch of the device, a 
stack of h-BN(46 nm)/Gr/h-BN(7 nm) on a Si/SiO2(300 nm) wafer used as global backgate. f. Photocurrent close to the resis¬ 
tance maximum at Vbg = —28 V shows a triangular photocurrent pattern, due to edge p-doping. The dashed lines in f indicate 
where the stack is underneath the gold, g. Photocurrent from charge puddles in encapsulated graphene on a metal gate close to 
the charge neutrality point (sketch in inset). The layers are AuPd(15 nm)/h-BN(42 nm)/Gr/h-BN(13 nm). The electrical contacts 
are on the left and right outside of this figure. In a-d,f,g the black solid lines indicate the graphene edge. 


ditions, gate voltages, and temperature, and evolve over 
timescales ranging from minutes to weeks. As an ex¬ 
ample, Figure 5a-d shows a progression of photocurrent 
maps, taken after annealing the device at 200 °C for 30 
minutes to temporarily remove charge density variations 
near the edges. Initially in Fig. 5a we see very small 
photocurrents indicating a flat carrier density landscape. 
After some time (^ hours), in the dark with only gate 
voltages smaller than 3 V applied, a small doping gradi¬ 
ent between the contacts builds up. This gradient leads 
to the stronger photocurrent shown in Fig. 5b. The lo¬ 
cal charge neutrality point, indicated by the maximum 
of photocurrent, is at the same position close to the edge 
of the device as further inside the bulk. After keeping 
the device for 3 hours in ambient conditions we can see 
a change of the local charge neutrality point at the edge 
of the graphene compared to the bulk in Fig. 5c. The 
edge is slightly more p-type compared to the bulk. Fi¬ 
nally we apply high gate voltages, of in this case 50 V 
for ^ 20 hours, to increase the edge doping. A strong 
p-doping at the edge and an n-doping in the bulk of the 
graphene is induced in Fig. 5d. This indicates that elec¬ 
tric field accelerates the speed and increases this type of 


edge doping. 

We exploit the observed edge doping to create a natural 
pn-junction along the edge of the device. For this we ap¬ 
ply a backgate voltage at which the edge of the graphene 
is p-type and the bulk n-type. We observe photocur¬ 
rent at the junction in Fig. 5f around the whole device, 
indicating that the edge doping is uniform around the 
graphene. The photocurrent decays gradually towards 
the midline between the electrodes as a result of how the 
triangular geometry modifies the ability of the contacts 
to capture photocurrents.We are able to temporarily 
reset the edge doping by annealing the device on a hot¬ 
plate at 200 °C for 30 minutes. 

While we have not been able to precisely identify the 
origin of the edge doping, we present here a technique to 
completely eliminate it. We place encapsulated graphene 
on top of a local conductive gate, such as a 15 nm AuPd 
alloy in the case of Fig. 5g. We find that edge doping is 
efficiently suppressed even after extended periods of time 
at ambient conditions and high gate voltages. Further¬ 
more, such devices lack the photodoping effect observed 
for devices where the h-BN is in contact with Si02^ (see 
Supplement, ID). We suspect that humidity that is able 
































6 


to penetrate between the boron nitride and the silicon 
oxide leading to trapped charges is responsible for the 
observed edge doping. 

In the device with a metal gate we find small features 
due to charge puddles on top of a slowly varying back¬ 
ground photocurrent, due to large scale carrier density 
inhomogeneities. The size of the features due to charge 
puddles determined by autocorrelation is ~ 800 nm. The 
long length scale of those features is either due to the 
longer cooling length of the encapsulated graphene com¬ 
pared to the graphene on Si 02 or due to larger charge 
puddle size in the encapsulated devices. Further work is 
required to clearly distinguish these effects. 

To conclude, we have demonstrated that scanning 
near-held photocurrent nanoscopy is a versatile tech¬ 
nique to characterize the electronic and opto-electronic 
and even previously inaccessible properties of relevant 
graphene devices. This technique is highly promising for 
spatially resolved quality control of regular graphene de¬ 
vices without the need for special device structures and 
can therefore be readily applied. 

METHODS 

Photocurrent model 

Photocurrent /pc in graphene as generated by the photo¬ 
thermoelectric effect and is described 



where R is the total resistance including graphene, contacts 
and circuitry, W the device width and x the current how di¬ 


rection. This is valid for rectangular graphene devices and 
special care needs to be taken for arbitrary shapes, such as in 
Fig. 5.^^ For the temperature prohle T(x) we consider that 
the heat spreads in two dimensions with heat sinking to lat¬ 
tice and substrate, producing a T(x) prohle described by a 
modihed Bessel function of the second kind, with a hnite tip 
size correction (see Supplement, II B). A 25 nm hnite tip size 
correction was used for all simulations. 


Measurement details 

The s-SNOM used was a NeaSNOM from Neaspec GmbH, 
equipped with a CO 2 laser operated at 10.6 pm, away from 
the phonon resonance of Si 02 , which can lead to strong sub¬ 
strate contributions to the photocurrent.^^ The probes were 
commercially-available metallized atomic force microscopy 
probes with an apex radius of ~ 25 nm. The tip height was 
modulated at a frequency of ~ 250 kHz with an amplitude 
of 60-80 nm. A Femto DLPCA-200 current pre-ampliher was 
used. 

Device fabrication 

The CVD graphene was transferred onto a self-assembled 
monolayer^^ on 285 nm of Si02 in order to stabilize the charge 
neutrality point. The contacts were dehned using optical 
lithography with Ti (5 nm)/Pd (35 nm). The graphene was 
transferred onto deposited contacts. 

The exfoliated graphene device was fabricated on 
a Si/SiO2(300 nm) wafer, used as backgate. The 
Cr(0.8 nm)/Au(80 nm) contacts were defined using electron 
beam lithography. 

The Si/SiO2(300 nm)/h-BN(46 nm)/Gr/h-BN(7 nm) and 
the Si/SiO2(300 nm)/AuPd(15 nm)/h-BN(42 nm)/Gr/h- 
BN(13 nm) stacks, were fabricated using the polymer-free van 
der Waals assembling technique. 


These authors contributed equally 

^ frank.koppens@icfo.eu 

^ Li, X. et al. Large-area synthesis of high-quality and uni¬ 
form graphene films on copper foils. Science 324, 1312- 
1314 (2009). 

^ Reina, A. et al. Large area, few-layer graphene films on 
arbitrary substrates by chemical vapor deposition. Nano 
Lett 9, 30-35 (2009). 

^ Bae, S. et al. Roll-to-roll production of 30-inch graphene 
films for transparent electrodes. Nature Nanotech. 5, 574- 
578 (2010). 

^ Bonaccorso, F. et al. Production and processing of 
graphene and 2d crystals. Mater. Today 15, 564-589 
( 2012 ). 

^ Ren, W. & Gheng, H.-M. The global growth of graphene. 
Nature Nanotech. 9, 726-730 (2014). 

® Lee, J.-H. et al. Wafer-Scale Growth of Single-Grystal 
Monolayer Graphene on Reusable Hydrogen-Terminated 
Germanium. Science 344, 286-289 (2014). 

^ Gao, L. et al. Face-to-face transfer of wafer-scale graphene 
films. Nature 505, 190-194 (2014). 

® de Heer, W. A. et al. Large area and structured epitaxial 


graphene produced by confinement controlled sublimation 
of silicon carbide. PNAS 108, 16900-16905 (2011). 

^ Ferrari, A. G. et al. Science and technology roadmap 
for graphene, related two-dimensional crystals, and hybrid 
systems. Nanoscale 7, 4598-4810 (2014). 

Akinwande, D., Petrone, N. & Hone, J. Two-dimensional 
flexible nanoelectronics. Nature Commun. 5, 5678 (2014). 
Koppens, F. H. L. et al. Photodetectors based on graphene, 
other two-dimensional materials and hybrid systems. Na¬ 
ture Nanotech. 9, 780-793 (2014). 

Yazyev, O. V. & Louie, S. G. Electronic transport in poly¬ 
crystalline graphene. Nature Mater. 9, 806-809 (2010). 
Yu, Q. et al. Gontrol and characterization of individual 
grains and grain boundaries in graphene grown by chemical 
vapour deposition. Nature Mater. 10, 443-449 (2011). 
Gummings, A. W. et al. Gharge transport in poly¬ 
crystalline graphene: challenges and opportunities. Adv. 
Mater. 26, 5079-5094 (2014). 

Yazyev, O. V. & Ghen, Y. P. Polycrystalline graphene 
and other two-dimensional materials. Nature Nanotech. 9, 
755-767 (2014). 

Martin, J. et al. Observation of electron-hole puddles in 



7 


graphene using a scanning single-electron transistor. Na¬ 
ture Phys. 4, 144-148 (2007). 

Chen, J.-H., Jang, C., Xiao, S., Ishigami, M. & Fuhrer, 
M. S. Intrinsic and extrinsic performance limits of 
graphene devices on Si02. Nature Nanoteeh. 3, 206-209 
(2008). 

Zhang, Y., Brar, V. W., Girit, C., Zettl, A. & Crommie, 
M. F. Origin of spatial charge inhomogeneity in graphene. 
Nature Phys. 5, 722-726 (2009). 

Decker, R. et al. Local electronic properties of graphene on 
a BN substrate via scanning tunneling microscopy. Nano 
Lett. 11, 2291-2295 (2011). 

Xue, J. et al. Scanning tunnelling microscopy and spec¬ 
troscopy of ultra-flat graphene on hexagonal boron nitride. 
Nature Mater. 10, 282-285 (2011). 

Burson, K. M. et al. Direct imaging of charged impurity 
density in common graphene substrates. Nano Lett. 13, 
3576-3580 (2013). 

Duong, D. L. et al. Probing graphene grain boundaries 
with optical microscopy. Nature 490, 235-239 (2012). 
Huang, P. Y. et al. Grains and grain boundaries in single¬ 
layer graphene atomic patchwork quilts. Nature 469, 389- 
392 (2011). 

Yasaei, P. et al. Bimodal Phonon Scattering in Graphene 
Grain Boundaries. Nano Lett. 15, 4532-4540 (2015). 
Deshpande, A., Bao, W., Miao, F., Lau, G. & LeRoy, B. 
Spatially resolved spectroscopy of monolayer graphene on 
Si02. Phys. Rev. B 79, 205411 (2009). 

Gibertini, M., Tomadin, A., Guinea, F., Katsnelson, M. I. 
& Polini, M. Electron-hole puddles in the absence of 
charged impurities. Phys. Rev. B 85, 201405 (2012). 

Gho, S. et al. Thermoelectric imaging of structural disorder 
in epitaxial graphene. Nature Mater. 12, 913-918 (2013). 
Fei, Z. et al. Electronic and plasmonic phenomena at 
graphene grain boundaries. Nature Nanoteeh. 8, 821-825 
(2013). 

29 Perrari, A. G. & Basko, D. M. Raman spectroscopy as 
a versatile tool for studying the properties of graphene. 
Nature Nanoteeh. 8, 235-246 (2013). 

Mueller, T., Xia, E., Ereitag, M., Tsang, J. & Avouris, 
P. Role of contacts in graphene transistors: A scanning 
photocurrent study. Phys. Rev. B 79, 245430 (2009). 
Mauser, N. & Hartschuh, A. Tip-enhanced near-held op¬ 
tical microscopy. Chem. Soe. Rev. 43, 1248-1262 (2014). 
Mauser, N. et al. Antenna-enhanced optoelectronic prob¬ 
ing of carbon nanotubes. Nano Lett. 14, 3773-3778 (2014). 
Grover, S., Dubey, S., Mathew, J. P. & Deshmukh, M. M. 
Limits on the bolometric response of graphene due to 
dicker noise. Appl. Phys. Lett. 106, 051113 (2015). 

Wang, L. et al. One-dimensional electrical contact to a 
two-dimensional material. Scienee 342, 614-617 (2013). 
Kretinin, A. V. et al. Electronic properties of graphene en¬ 
capsulated with different two-dimensional atomic crystals. 
Nano Lett. 14, 3270-3276 (2014). 

Ju, L. et al. Photoinduced doping in heterostructures of 
graphene and boron nitride. Nature Nanoteeh. 9, 348-352 
(2014). 

Eei, Z. et al. Infrared nanoscopy of dirac plasmons at the 
graphene-Si02 interface. Nano Lett. 11, 4701-4705 (2011). 
Keilmann, E. & Hillenbrand, R. Near-held microscopy by 
elastic light scattering from a tip. Phil. Trans. R. Soe. A 
362, 787-805 (2004). 

Lemme, M. G. et al. Gate-Activated Photoresponse in a 
Graphene p-n Junction. Nano Lett. 11, 4134-4137 (2011). 


Gabor, N. M. et al. Hot carrier-assisted intrinsic photore¬ 
sponse in graphene. Seienee 334, 648-652 (2011). 

Song, J. G. W., Rudner, M. S., Marcus, G. M. & Levitov, 
L. S. Hot carrier transport and photocurrent response in 
graphene. Nano Lett. 11, 4688-4692 (2011). 

Tielrooij, K. J. et al. Hot-carrier photocurrent ehects at 
graphene-metal interfaces. J. Phys. Condens. Matter 27, 
164207 (2015). 

Badioli, M. et al. Phonon-Mediated Mid-Infrared Photore¬ 
sponse of Graphene. Nano Lett. 14, 6374-6381 (2014). 
Song, J. G. W. & Levitov, L. S. Shockley-Ramo theorem 
and long-range photocurrent response in gapless materials. 
Phys. Rev. B 90, 075415 (2014). 

Zuev, Y., Ghang, W. & Kim, P. Thermoelectric and Mag¬ 
netothermoelectric Transport Measurements of Graphene. 
Phys. Rev. Lett. 102, 096807 (2009). 

Tsen, A. W. et al. Tailoring electrical transport across 
grain boundaries in poly crystalline graphene. Seienee 336, 
1143-1146 (2012). 

Van Tuan, D. et al. Scaling properties of charge trans¬ 
port in polycrystalline graphene. Nano Lett. 13, 1730-1735 
(2013). 

Wei, R, Bao, W., Pu, Y., Lau, G. N. V Shi, J. Anomalous 
thermoelectric transport of dirac particles in graphene. 
Phys. Rev. Lett. 102, 166808 (2009). 

Lee, E. J. H., Balasubramanian, K., Weitz, R. T., 
Burghard, M. & Kern, K. Gontact and edge ehects in 
graphene devices. Nature Nanoteeh. 3, 486-490 (2008). 
Shalom, M. B. et al. Proximity superconductivity in bal¬ 
listic graphene, from Eabry-Perot oscillations to random 
Andreev states in magnetic held. arXiv 1504.03286 (2015). 
Allen, M. T. et al. Spatially resolved edge currents 
and guided-wave electronic states in graphene. arXiv 
1504.07630 (2015). 

Ghen, S. Y., Ho, P. H., Shine, R. J., Ghen, G. W. X Wang, 
W. H. Transport/magnetotransport of high-performance 
graphene transistors on organic molecule-functionalized 
substrates. Nano Lett. 12, 964-969 (2012). 


ACKNOWLEDGEMENTS 

It is a great pleasure to thank Stijn Goossens, Klaas-Jan 
Tielrooij and Misha Eogler for many fruitful discussions and 
Eabien Vialla for the assistance in rendering the graphics in 
Eig. la. Open source software was used (www.matplotlib.org, 
www.python.org,www.povray.org). E.H.L.K. acknowledges 
support by the Eundacio Gellex Barcelona, the ERG Ga- 
reer integration grant 294056 (GRANOP), the ERG starting 
grant 307806 (GarbonLight). E.H.L.K. and R.H. acknowledge 
support by the E.G. under Graphene Elagship (contract no. 
GNEGT-IGT-604391). D.J. acknowledges support from the 
Ramon y Gajal fellowship program. Y.G. and J.H. acknowl¬ 
edge support from the US Office of Naval Research N00014- 
13-1-0662. We acknowledge financial support from the Span¬ 
ish Ministry of Economy and Gompetitiveness and ’’Eondo 
Europeo de Desarrollo Regional” through Grant TEG2013- 
46168-R. Q.M. and P.J.H. have been supported by AEOSR 
Grant No. EA9550-11-1-0225 and the Packard Eellowship pro¬ 
gram. This work made use of the Materials Research Science 
and Engineering Genter Shared Experimental Eacilities sup- 


ported by the National Science Foundation (NSF) (Grant No. 
DMR-0819762) and of Harvard’s Center for Nanoscale Sys¬ 
tems, supported by the NSF (Grant No. ECS-0335765). S.R. 
acknowledges the Spanish Ministry of Economy and Com¬ 
petitiveness for funding (MAT2012-33911), the Secretaria de 
Universidades e Investigacion del Departamento de Econo- 
mia y Conocimiento de la Generalidad de Catalunya and the 
Severo Ochoa Program (MINECO SEV-2013-0295). J.E.B.V. 
acknowledges support from SECITI (Mexico, D.E.). 

AUTHOR CONTRIBUTIONS 

A.W. performed the experiments, analysed the data and 
wrote the manuscript. P.A.-G. helped with measurements, 
interpretation and discussion of the results. M.B.L. helped 
with data analysis, measurements, interpretation, discussion 


of the results and manuscript writing. Y.G. fabricated the 
h-BN/graphene/h-BN devices. G.N. and D.J. fabricated 
the CVD graphene devices. Q.M. fabricated the exfoliated 
graphene devices. K.W. and T.T. synthesized the h-BN. 
J.E.B.V., A.W.C. and S.R. performed the simulations of the 
Seebeck coefficient at grain boundaries. R.H. and E.H.L.K. 
supervised the work, discussed the results and co-wrote the 
manuscript. All authors contributed to the scientific discus¬ 
sion and manuscript revisions. 

COMPETING FINANCIAL INTERESTS 

R.H. is co-founder of Neaspec GmbH, a company produc¬ 
ing scattering-type scanning near-held optical microscope sys¬ 
tems such as the ones used in this study. All other authors 
declare no competing hnancial interests. 



Supplementary Material: Near-field photocurrent nanoscopy on bare 

and encapsulated graphene 


1 


I. MEASUREMENTS 


In this section we show measurements that further support the proof that the antisymmetric near-held photocurrent 
pattern of Fig. 1 of the main text indeed stems from a grain boundary. 


A. Topography, near-held optical scattering, and near-held photocurrent near a grain boundary 


In Fig. Ic and d of the main text, which are reprinted here as Fig. Sla and c, we show the topography and 
the near-held photocurrent at a grain boundary for a backgate voltage of 0 V. To provide further evidence that the 
photocurrent stems from a grain boundary, we also show the simultaneously acquired near-held optical scattering data 
in Fig. Sib. The near-held optical scattering data is acquired using a conventional scattering-type scanning near-held 
optical microscope (s-SNOM).^^ It is a measure of the light that is scattered out of the graphene that interacts with 
the tip and is consequently detected in the far-held. The near-held optical scattering data show the typical double 
fringe due to plasmon rehections at the grain boundary. 


a b c 



Figure SI | Comparison between topography, near-field optical scattering and near-field photocurrent, a, Topography of 
etched CVD graphene does not show grain boundaries but only wrinkles and other inhomogeneities due to the transfer process, 
b, Near-field optical scattering shows the characteristic plasmonic double fringes around a grain boundary in CVD graphene at 
Vbg = 0 V because the carrier density is ns = 3.7 x 10^^ cnn“^ and plasmons are supported.c, Near-field photocurrent clearly 
shows the grain boundary and a sign change around it. The dashed line in b,c indicates where the backgate dependent measure¬ 
ments of the near-field optical scattering and near-field photocurrent in Figs. 2c and d respectively in the main text where taken. 
All measurements of the near-field optical scattering ^opt presented here were obtained from the third harmonic interferometric 
pseudo-heterodyne signal,measured with a cryogenic HgCdTe detector. For simplicity Fig. Sib only shows |^opt|- 


Even though conventional mid-infrared s-SNOM also offers the ability to detect grain boundaries in the near-held 
optical scattering,there are some distinct advantages to using the near-held photocurrent. In order to observe 
plasmons due to grain boundaries, the graphene needs to be highly doped, as plasmons in graphene only propagate 
at at elevated carrier densities.For small carrier densities plasmons are heavily damped and no double fringe 
is visible, as observed by a decreasing visibility of the double fringe near-held optical scattering close to the charge 
neutrality point of the graphene (see Fig. S2b). 

In contrast, the near-held photocurrent shows a clear signature of the grain boundary even for small carrier densities 
(Fig. S2c), which enables us to extract more information, as discussed in the main text. This demonstrates further 
that near-held photocurrent is a useful tool to provide insight into the local properties of graphene. The near-held 
photocurrent technique also facilitates measurements compared to measuring the near-held optical scattering as no 
weak outscattered light needs to be collected and no interferometric measurement is required. 

When comparing the topography in Fig. Sla with the near-held optical scattering in Fig. Sib and the near-held 
photocurrent in Fig. Sic it becomes evident that the near-held optical scattering and the near-held photocurrent 
contain information that is not visible in the topography. 
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Figure S2 | Gate dependent measurement of a grain boundary, a, Backgate dependence of the resistance of the device mea¬ 
sured simultaneously to the near-field optical scattering and near-field photocurrent, b, Near-field optical scattering shows the plas- 
mon reflections due to the grain boundary and their dependence on the carrier density of the graphene, c, Near-field photocurrent 
shows no sign change even for highly doped graphene. The dashed dotted curves in b and c show the theoretical fringe spacing for 
a phase shift due to the reflection of —3/47r.^^ 



B. Gate dependence of a grain boundary 

As explained in the main text and shown in Fig. 2c, the Seebeck coefficient at or very near the grain boundary is 
smaller in magnitude than the Seebeck coefficient of the surrounding pristine graphene for all the carrier densities 
measured. To further support this statement, we show the near-held optical scattering and near-held photocurrent 
for an extended carrier density range in Fig. S2. These data show that there is no additional sign change for higher 
carrier densities indicating that the Seebeck coefficient at the grain boundary 5 'gb is smaller in magnitude than the 
Seebeck coefficient of pristine graphene Sq for the measurable range of carrier densities. The measurements where 
done on the same device and grain boundary as the ones shown in the main text. We note that in the case of Fig. S2c 
there is no sign change as a function of gate voltage visible as the charge neutrality point of the device was not reached 
due to high intrinsic doping. 

When comparing the position of the photocurrent extrema in Fig. S2c with the expected plasmonic fringe spacing 
as indicated by the dashed dotted curve it becomes evident indeed also the position of the photocurrent extrema is 
changing with carrier density. This could be due to an increased absorption due to the excitation of plasmons in the 
graphene and is subject of further study. 


C. Photodoping and puddles of encapsulated graphene 



Figure S3 | Encapsulated graphene device on oxidized silicon wafer, a, Near-field photocurrent of the graphene triangle de¬ 
vice at -90 V backgate voltage in the dark, b, Photodoping showing up in the near-field photocurrent after illuminating the device 
for several minutes with white LED light at —90 V. The dashed lines in a and b indicate where the stack is underneath the gold 
and the solid lines the outer edge of the stack and gold respectively. 


It is known that graphene on h-BN on top of an oxidized silicon wafer shows photodoping when it is illuminated 
with visible light.This effect is clearly visible in our near-field photocurrent measurements after illumination with 
visible light. In Fig. S3a we see almost no near-held photocurrent at a backgate voltage of —90 V. In Fig. S3b we 
show the near-held photocurrent pattern after the sample was illuminated with a white LED light source for several 
minutes at a backgate voltage of —90 V. Clearly, near-held photocurrent is visible all the way throughout the device. 
We attribute this to the screening of the backgate by photoexcited defects in the h-BN or at the h-BN/Si02 interface. 
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which can effectively neutralize the graphene.Thus we see charge puddles which show a very similar near-field 
photocurrent pattern compared to the near-field photocurrent from charge puddles observed at the charge neutrality 
point of exfoliated graphene on Si02 shown in Fig. 3 and Fig. 4 of the main text. In the case of the encapsulated 
device the charge puddles are induced by the photoexcited charged defects in the h-BN. In order to reset the photo 
doping and be immune to it we use positive gate voltages. 


D. Encapsulated graphene with a local gate 
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Figure S4 | Encapsulated graphene with independently tunable gates, a, 6-fold pattern showing the thermoelectric origin of 
the near-field photocurrent at the junction between the two gates.b, Near-field photocurrent from a local gated encapsulated 
graphene device with both sides of the local gate being tuned to a voltage of Vbg = 0.4 V close to the charge neutrality point. 
The layers are AuPd(15 nm)/h-BN(42 nm)/Gr/h-BN(13 nm). The dashed line indicates the position of the 50 nm gap between 
the two local gates. 


We also studied near-field photocurrent from an encapsulated device which was put onto conductive PdAu alloy 
gates with a 50 nm gap in between them in order to individually tune the carrier density in the graphene above the 
two gates. In Fig. S4a we show a six-fold pattern typical for a photothermoelectric effect^^ in graphene measured at 
the junction between the two gates, which is indicated by the dashed line in Fig. S4. 

In Fig. S4b the current in the device is shown for both sides being tuned to a gate voltage of 0.4 V, close to the 
charge neutrality point. The magnitude of the near-field photocurrent is two orders of magnitude smaller than the 
magnitude of near-field photocurrent from charge puddles in graphene on Si02. This indicates that the charge density 
inhomogeneity of the charge puddles in encapsulated graphene on a local gate is much lower than for graphene on 
Si02. Furthermore the length scale of the charge puddles indicates that the charge puddle size is much larger than 
for the case of graphene on Si 02 ^^^’^^^ or that the cooling length in encapsulated graphene devices is much longer 
than for bare graphene, as expected due to the increased carrier mobility in encapsulated graphene. 

We remark the we do not see any photodoping for the encapsulated graphene on top of a conductive gate structure, 
even after extended periods of exposure of many minutes to the same white light LED as was used to induce photodop¬ 
ing in Fig. S3b. This can be explained by the extraction of the photoexcited charged defects by the conductive gate, 
which is in direct contact with h-BN. 


II. MODELLING 


In this section we give a more detailed overview of the photothermoelectric photocurrent model used to describe 
the measurements and models in the main text. 


A. Photothermoelectric model 


Near-held photocurrent Jpc in graphene is governed by the photothermoelectric effect, which can be calculated 

Py.S9,S12,S13 


Ipcix,y) = ^ j ^^-^j^Six,y)dxdy 


(si) 
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where R is the total resistance, consisting of graphene, contact and circuit resistance, W the width of the device, 
dTjdx is the gradient of the temperature T in current flow direction x and S is the Seebeck coefficient. In order to 
simulate the near-held photocurrent in arbitrary geometries special care has to be taken. 

For calculating /pc at each position one has to do a convolution of the spatial Seebeck coefficient prohle with the 
temperature gradient. This is numerically very expensive and the simulations can take a long time. In order to 
more efficiently simulate a near-held photocurrent map for an arbitrary spatial Seebeck prohle we use the convolution 
theorem. This allows us to just multiply the Fourier transform of the spatial Seebeck coefficient and temperature 
gradient, multiply them and inverse Fourier transform them. This is computationally much faster and the simulation 
time is greatly reduced. This method was employed for all simulations shown in the main text as well as the 
Supplement. 


B. Heat spreading in 2D 

In order to describe the near-held photocurrent in graphene with a model it is of great importance to correctly 
describe the heat prohle within the graphene as the heat prohle in combination with the Seebeck prohle lead to the 
near-held photocurrent pattern. 

Graphene on top of a substrate is a two dimensional material with the substrate acting as a heat sink, as for typical 
substrate materials such as Si02 the thermal conductivity is much larger than for air. The two dimensional heat 
equation in steady state with an additional term for heat sinking can be written as: 

0 = + P - g{T -T^), (s2) 

where T is the electron temperature, Tg is the constant heat sink temperature, n the thermal conductivity in plane, 
g the thermal conductivity out of plane to the heat sinking substrate and P the power density of the heat source. 



Figure S5 | Steady state temperature profile created in infinite graphene with a thermal length of 200 nm and a finite tip size 
correction of 25 nm. 

We apply this equation by assuming that the electron heat cools to the lattice, but then that the lattice can pass 
on any heat much more easily to substrate. In that way the temperature of the lattice does not change significantly. 
This is convenient because the lattice itself has a different thermal length which would complicate matters. 

If we take the two-dimensional case and a point source for P, 


p = PtotaiS{x)S{y), (s3) 

then the solution of the two dimensional heat eq. s2 turns out to be: 

T,pot =Toi^o (s4) 

V Gool / 


where Kq is the modified Bessel function of the second kind, r = the cooling length /cool = and 

the maximum temperature rise Tq = Ptotai/(27r/^). 

If we also include /tip to approximate the effect of the finite tip size we end up with a heat spot of the following 
form: 


Fspot {x' -x,y -y)= TqKo 


I^V ^cool 


(s5) 


A typical temperature profile as calculated by eq. s5 with a finite tip size approximation /tip = 25 nm, the actual 
radius of the tip, and a cooling length of /cool = 200 nm is presented in Fig. S5. 
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C. Grain boundary model 

Grain boundaries can be modelled as having a finite width with a Gaussian profile^^ and their Seebeck coefficient 
is smaller in magnitude than the one of the surrounding pristine graphene, in accordance with the results of the main 
text. The near-held photocurrent prohle is then calculated by performing a two dimensional convolution between the 
temperature prohle dehned by eq. (s5) and the Seebeck prohle. The results of this convolution are shown in Fig. S6. 
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Figure S6 | Simulated near-field photocurrent map from a grain boundary. A map of the simulated near-field photocurrent 
for each position of the graphene device. The contacts are on the left and right outside of the region shown in this figure. The 
cooling length used was 140 nm and the tip size 25 nm. 


D. Seebeck coefficient of polycrystalline graphene 

We use an order-N Kubo-Greenwood wavepacket approach to calculate the conductivity of polycrystalline samples 
with diherent average grain sizes,and we use square samples to convert this conductivity to conductance (G). 
Additionally, we calculate the conductance of pristine graphene using a Landauer approach. The Seebeck coefficient 
is calculated as the ratio of the first- and zero-order Onsager coefficients. 


S{fi,T) = - 


1 C,oiE-lAG{E)(^-§i'jdE 
^ fZG(E)(-§i)‘lE 


(s6) 


where / is the Fermi-Dirac distribution, fi the chemical potential and we set the temperature T to 300 K. 

For polycrystalline samples with different average grain size (13, 18, 21 and 25 nm) we find that the Seebeck 
coefficient is significantly reduced compared to the clean case due to scattering at the grain boundaries (Figure S7). 
Furthermore, the Seebeck coefficient is independent of grain size, resulting from the linear scaling of charge transport 
in polycrystalline graphene^^^’^^^. The impact of the grain boundaries is reduced for larger chemical potentials 
(> 0.3 eV), where the Seebeck coefficient for polycrystalline and pristine graphene is similar. 



Figure S7 | Simulated Seebeck coefficient for different polycrystalline samples {T — 300 K), for grain sizes (13, 18, 21 and 25 
nm). 

The simulation of the backgate dependent near-held photocurrent in Fig. 2e of the main text was performed by 
using a Gaussian by using a Gaussian spatial distribution of the Seebeck coefficient with a full width at half maximum 
of 20 nm. The difference in Seebeck coefficient between Sq of the pristine graphene and 5 'gb at the grain boundary 
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was calculated for each measured backgate voltage. For 5 'gb, the simulation corresponding to 25 nm polycrystalline 
graphene was used, but the results are essentially the same for different grain sizes. The spatial near-held photocurrent 
prohle for each backgate voltage was calculated according to the procedure given in IIC and then normalized by the 
resistance R for each gate voltage according to eq. (si). 


E. Charge puddle model 


In order to model near-held photocurrent from charge puddles, we hrst need to hnd an accurate model of the charge 
puddle distribution. For this we use a random spatial Seebeck prohle generated by a spatial prohle of white noise and 
smoothing the noise to create charge puddles with an average approximate size of 20 nm. This size was extracted 
from previous experimental studies of charge puddles on 

We calculate the near-held photocurrent map from the charge puddles. To this end, we again convolve the spatial 
Seebeck prohle with the spatial heat gradient prohle. A comparison between the spatial Seebeck prohle and the 
near-held photocurrent map is shown in Fig. S8. Here the source and drain contacts are outside of the display on the 
left and right respectively. It is obvious that for positions with a high gradient in Seebeck coefficient the near-held 
photocurrent is strongest, as expected from eq. si. 
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Figure S8 | Simulation of a near-field phtocurrent map induced by a random charge puddle distribution, a, Randonn dis¬ 
tribution of Seebeck coefficients due to a randonn distribution of charge puddles with an average size of ~ 20 nm. b, Near-field 
photocurrent generated by the charge puddles in a with a cooling length of the charge carriers in the graphene of 200 nm and a tip 
size of 25 nm. 


F. Model of a pnp-junction 


Here, we model the near-held photocurrent map of a pnp-junction, where the p-type and n-type regions are extended 
to a size larger than the heat spot size. The transition length scale of the regions is smaller than the heat spot size. In 
this case the near-held photocurrent at the pn- and np-junction respectively has a single sign as observed for example 
in the triangular near-held photocurrent pattern due to edge doping in Fig. 5f in the main text. In the simulations 
shown in Fig. S9a and b the electrical contacts are on the left and right outside of the shown region. 
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Figure S9 | Simulation of a near-field photocurrent map from a pnp-junction. a, Spatial Seebeck coefficient profile for two 
p- and one n-type regions. The doping in the two regimes is assumed to be of the same magnitude and opposite sign, b, Near-field 
photocurrent generated by the pn- and np-junction with a cooling length of 200 nm and a tip size of 25 nm. 
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